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We report a study of overdamped Josephson junction arrays with the geometry of periodically 
repeated Sierpinski gaskets. These model superconductors share essential geometrical features with 
truly random (percolative) systems. When exposed to a perpendicular magnetic field B, their 
euclidian or fractal behavior depends on the relation between the intervortex distance (imposed by 
B) and the size of a constituent gasket, and was explored with high-resolution measurements of the 
sample magnetoinductance L(B). In terms of the frustration parameter / expressing (in units of 
the superconducting flux quantum) the magnetic flux threading an elementary triangular cell of a 
gasket, the crossover between the two regimes occurs at f C N = 1/(2 x 4^), where N is the gasket 
order. In the fractal regime (/ > / c jv) a sequence of equally spaced structures corresponding to 
the set of states with unit cells not larger than a single gasket is observed at multiples of f C N, as 
predicted by theory. The fine structure of L(f) radically changes in the euclidian regime (/ < / c jv), 
where it is determined by the commensurability of the vortex lattice with the effective potential 
created by the array. Anomalies observed in both the periodicity and the symmetry of L(f) are 
attributed to the effect of a hidden incommensurability, which arises from the deformation of the 
magnetic field distribution caused by the asymmetric diamagnetic response of the superconducting 
islands forming the arrays. 

PACS numbers: 74.80.-g, 74.50.+r, 74.60. Ge, 74.25.Nf 



I. INTRODUCTION 

A wide variety of disordered materials, including su- 
perconductors, is known to exhibit geometrical inhomo- 
geneities over a broad range of length scales. The prop- 
erties of such systems can be conveniently described in 
terms of percolationu, the simplest idea to understand 
randomness. Percolation can be regarded as a geomet- 
rical phase transition taking place at a "critical concen- 
tration" p c separating a phase of finite clusters (p < p c ) 
from a phase where an infinite cluster is present (p > p c ). 
Like other critical phenomena, it is characterized by a 
correlation length £ p (p) which diverges at the percola- 
tion threshold p c . Right at p c , a system with percolative 
disorder exhibits a natural self-similar structure at all 
length scales and can therefore be modeled by a family 
of scale-invariant fractal lattices, such_as the Sierpinski 
gasket (SG), which has been suggested^ to mimic the es- 
sential geometrical features of the percolating cluster's 
backbone. In the critical region above and below p c , 
where £ p is finite, the nature of the geometry depends on 
the length scale I at which one is probing the system: if 
I < £ p , its structure is fractal, whereas it can be regarded 
as homogeneous with conventional euclidian features for 

Allowing an accurate control of both the nature and 
the amount of disorder and exhibiting properties quite 
sensitive to dimensionality, Josephson junction arrays 
and wire networks prepared with modern micro- and 
nanofabrication techniques provide ideal model systems 
to investigate the dimensional crossover from the euclid- 



ian to the fractal regime. The .first step in this direc- 
tion was made by Gordon et alu, who investigated the 
superconducting-to-normal phase boundary T C (B) of alu- 
minum wire networks formed by periodically repeated 
SGs and of analogous networks with percolative geometry 
exposed to a perpendicular magnetic field B. In those ex- 
periments the magnetic length, 1(B) w (tfio/B) 1 / 2 , which 
is a measure of the typical nearest-neighbor distance be- 
tween the vortices present in the system, was shown to be 
the relevant length scale to explore the euclidian-fractal 
(EF) crossover (</> is the superconducting flux quantum). 
While the scaling behavior of the phase boundary of the 
SG networks was found to exhibit a crossover from the 
euclidian to the fractal regime consistent with theoreti- 
cal predictions!] based on extensionsOu of theJGIinzburg- 
Landau theory and allowing for comparison^ with the 
anomalous diffusion exponent, no EF crossover was ob- 
served in percolative networks which, surprisingly, were 
found to behave like a homogeneous system at all length 
scales 1(B) probed in the experiment. 

More recently, compelling evidence for the EF 
crossover in a percolati'va system emerged from ac con- 
ductance measurementsu performed in zero magnetic 
field on (unfrustrated) site-diluted Josephson junction ar- 
rays with site occupation probabilities p very close to p c . 
In these experiments the crossover was controlled by the 
driving angular frequency u>, which determines the ratio 
of the impedances associated with the two types of links 
forming the random network. According to Efros and 
ShklovskiiEa, the increase of this ratio with decreasing uj 
also leads, if p is sufficiently close to p c , to a crossover 
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FIG. 1. (a) Scanning electron micrograph showing a por- 
tion of a triangular array of periodically repeated second-order 
Sierpinski gaskets of proximity-effect coupled Pb/Cu/Pb 
Josephson junctions. The length of the elementary links of 
the gaskets is 8/xm. Notice the "truncated-star" shape of the 
superconducting Pb islands (with the exception of those cen- 
tered at the common vertices of three constituent gaskets), 
(b) The rhombohedral unit cell of a periodic array of sec- 
ond-order gaskets. 

from the fractal to the euclidian regime, as confirmed by 
the experiments of Ref. ||. 

In this article we report the results of experimental 
and theoretical studies of a model superconductor shar- 
ing essential geometrical features with a percolative sys- 
tem. The samples we have investigated are arrays of 
proximity-effect coupled SNS junctions (where S stands 
for superconductor and N for normal metal) consisting 
of SGs connected to each other at the vertices in such a 
way as to form a regular triangular lattice. As can be 
seen from Fig. 1(a), where a part of an array of second- 
order gaskets is shown, in these systems the linear size 
ijv = 2^0 of an individual gasket of order N can be 
regarded as playing the role of £ p (a is the length of an 
elementary link of the gasket). As shown in detail in 
this work, a remarkable feature of the SG arrays is that 
in these systems the EF crossover is clearly manifest, in 
contrast to truly percolative systems, where its signatures 
are elusive. 

The quantity at the heart of the present study is the 
magnetoinductancc L(B) of the SG arrays, extracted 
from measurements of their ac impedance. Its inter- 
est resides in the observation that, being inversely pro- 
portional to the areal superfluid density, it provides a 
tool to appreciate how the degree of superconducting 
phase coherence in the system changes with B or, equiva- 
lently, with the level of frustration imposed by B. Previ- 
ous impedance measurements performed on weakly frus- 
trated arrays similar to those studied in this work fo- 
cused merely on the fractal regime and demonstrated, in 
particular, the unusual scaling properties£3 of the vor- 



tex energyO as well as the asymptotic (B — > 0) scaling 
behavior of the field-induced correction to the array in- 
ductance, resulting from the hierarchical structure of the 
gasketalj. 

Although some preliminary evidence for a dimensional 
crossover was already reported in Ref. [l2], the phe- 
nomenon was not exhaustively investigated. In the 
present work we rely on high-resolution studies of the 
complex fine structure of L(B) to explore in detail both 
the euclidian and the fractal regimes of the SG ar- 
rays. Reflecting flux quantization phenomena occurring 
in loops with a hierarchical distribution of sizes up to the 
gasket size Ln, the fine structure provides a unique tool 
to reveal how the geometrical properties of the system 
change as the magnetic length 1(B) is swept through Ln- 
We show that, in terms of the frustration parameter / ex- 
pressing the magnetic flux threading an elementary trian- 
gular cell of a gasket in units of (j>o, the EF crossover oc- 
curs at fcN — 1/(2 x 4 7V ). In the fractal regime (/ > f cN ) 
the most relevant contributions to the fine structure of 
£(/) are shown to arise from a particular set of ground 
states defined by / = M/ C jy, where M is an integer. 
Corresponding to vortex configurations where the vor- 
tex lattice is strongly pinned by the hierarchical poten- 
tial landscape created by the gaskets, these states are 
particularly robust against thermal fluctuations and are 
therefore quite prominent in the fine structure of L(f). 

A very interesting aspect emerged from the study of 
the magnetoinductance in the fractal regime. The analy- 
sis of the data revealed anomalous features (specifically, 
the suppression of the periodicity corresponding to a shift 
of / by 1 and of the symmetry with respect to / = 1/2) 
inconsistent with theoretical predictions based on the de- 
scription of the system in terms of a uniformly frustrated 
XY model. We suggest that, because of the asymmetric 
shape of the superconducting islands forming the junc- 
tion pattern of a gasket [see Fig. 1(a)], the screening cur- 
rents flowing in the islands create a distortion of the mag- 
netic field distribution in the array such that the fluxes 
threading the various loops slightly deviate from being 
proportional to their areas. This introduces an effective 
incommensurability (which we call "hidden" to distin- 
guish it from the "geometric" one studied earlierLj in 
wire networks with incommensurate cells) and perturbs 
the self-similarity of the gaskets. As a result, the system 
is no longer uniformly frustrated. We demonstrate that 
the anomalous features mentioned above can be quanti- 
tatively accounted for by a simple model, in which the 
area of the different plaquettes of a gasket is changed 
according to an appropriate deformation scheme. 

Compelling evidence for the existence of the euclidian 
regime is provided by the array's magnetoinductance for 
/ < fcN- Besides the absence of the power-law scal- 
ing behavior characteristic of the fractal regime, L(f) 
contains structures which reflect the presence of ground 
states corresponding to vortex configurations commen- 
surate with the underlying triangular lattice formed by 
the largest triangular loops of the array (below f C N , it is 
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energetically unfavorable for vortices to penetrate loops 
of smaller size) , thereby allowing an unambiguous iden- 
tification of the euclidian regime. 

In general, these vortex configurations have the same 
structure as in the XY model on a honeycomb lattice, 
but correspond to values of / reduced by a factor of 
l//cJV = 2x4 N . However, the case / = (1/2) f C N requires 
special attention. In a honeycomb lattice the ground 
state of the corresponding (fully frustrated) XY model 
is characterized by an accidental degeneracy, which (to 
lowest ordep)i survives even in the presence of thermal 
fiuctuationaU. Owing to the more complex structure of 
the system, however, this peculiar degeneracy is removed 
in our SG arrays, allowing to identify the vortex config- 
uration in the ordered phase at / = (1/2)/ c jv- 

The paper is organized as follows. Experimental de- 
tails are given in Sec. II. Sec. Ill is devoted to the 
fractal regim e ( f > / c jv). Relying on the methods devel- 
oped in Ref. [jJL in Sec. Ill A we present the calculation 
of the magnetoinductance of a frustrated SG array for 
the particular set of frustrations / = Mf c ^ correspond- 
ing to the sequence of the most stable states, which are 
characterized by a relatively compact structure (with a 
unit cell not larger than a single gasket). Experimental 
data for the magnetoinductance in the same regime are 
presented and discussed in Sec. Ill B, where we show, 
in particular, how the anomalous features of L(f), re- 
vealing the presence of hidden incommensurability, can 
be accounted for by a simple model, in which the areal 
changes of different plaquettes (related to the redistribu- 
tion of the magnetic field) are determined by only one 
adjustable parameter. In Sec. IV we provide experimen- 
tal evidence for the existence of the EF crossover and in 
Sec.V the fine structure of L(f) in the euclidian regime 
below f c ]y is shown to be consistent with the existence 
of vortex states commensurate with the periodic lattice 
formed by the largest triangular cells of the arrays. A 
few concluding remarks are given in Sec. VI. 

II. EXPERIMENTAL ASPECTS 

The samples studied in this work consist of second- 
(N = 2) and fourth-order (N — 4) gaskets sitting on the 
sites of, respectively, a313x313 and a 78 x 78 triangular 
lattice and connected to each other at the vertices [see 
Fig. 1(a)]. Each gasket contains, respectively, 3 2+1 = 
27 and 3 4+1 = 243 SNS Josephson junctions consisting 
of superconducting Pb islands proximity-effect coupled 
to each other by an underlying normal Cu layer. The 
geometrical and physical parameters of the junctions are 
almost identical to those of the array studied in Ref. |lj. 
Most of the data presented below have been obtained 
in experiments performed on the array of second-order 
gaskets. 

The sheet magnetoinductance L(f) was inferred from 
measurements of the array's linear sheet impedance Z = 



R+iluL performed with a very sensitive-SQUID-operated 
two-coil mutual inductance techniquet^l at driving fre- 
quencies typically in the range 0.1-1.0 kHz. With this 
method we were able to resolve inductance changes of 
the order of 10 pH in swept-frustration measurements, 
in which / could be tuned with a precision better than 
10~ 3 . The experimental data are presented and analyzed 
in terms of L _1 (/), the quantity measuring the degree of 
superconducting phase coherence in the samples. When 
needed, the resistive component R(f), related to dissipa- 
tive vortex motion, is also shown for completeness. Addi- 
tional details concerning the samples and the measuring 
technique can be found in Ref. [l^. 

In the following, temperatures are expressed in 
terms of the reduced temperature relevant for the sta- 
tistical mechanics of the system, r = fcgT/ J(T), 
where J(T) is the temperature-dependent Josephson 
coupling energy. At temperatures well below the 
zero-field critical temperature r c jv, J(T) was deduced 
from measurements of the "bare" sheet kinetic induc- 
tance L(T) =U /2tt) 2 (5/3) n /V3J(T) of the unfrus- 
trated samplesO. Extrapolation to higher temperatures 
was then achieved by fitting the low-temperature data to 
theoretical expressionsE3 for J(T). 

Because of their two-dimensional (2D) nature at length 
scales larger than Ln, both samples are expected to 
exhibit, at zero frustration (/ = 0) and in the limit 
uj — * 0, a-J3crczinskii-Kosterlitz-Thouless (BKT) phase 
transitiorO. A sharp depression of L _1 (0), which can 
be associated with the BKT transition, has been indeed 
observed at, respectively, r c2 « 0.57 and t C 4 Ri_£l.23, in 
good agreement with the theoretical predictiorcHl t c n — 
{3/5) N Tco, where r c0 ~ 1.5 is the reduced temperature 
of the BKT transition of a regular triangular Josephson 
junction arraylij with the same J(T). 

III. THE FRACTAL REGIME 

A. The ground states of a regular array of Sierpinski 
gaskets and their sheet inductance 

We start by recalling that, within the framework of an 
approximation ignoring thermal fluctuations, a Joseph- 
son junction array behaves, with respect to an external 
(dc) current source, like a netw ork of inductors {Lij}, 
whose inductances are giver£3ll3 by: 

where 9ij is the gauge-invariant phase difference across 
the link ij. As required by fluxoid quantization, the sum 
of 6ij around a lattice cell is equal to 2ir(fS — to), where 
/ is the frustration parameter expressing the magnetic 
flux (in units of </>o) threading an elementary triangular 
cell of a gasket [/ = Ba 2 v / 3/(4</) )], S the area of the cell 
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(expressed in units of the area of an elementary trian- 
gular loop) and m the number of flux quanta (vortices) 
penetrating the cell under consideration. 

In writing Eq. (1) we have assumed that the proximity- 
effect coupled SNS junctions forming the arrays studied 
in this work have a sinusoidal current-phase relation at 
the temperatures of interested. It clearly follows from Eq. 
(1) that, even if all the junctions are identical, their ef- 
fective inductances in a frustrated system may differ sub- 
stantially from each other on account of the nonuniform 
distribution of {Oij}. Since the array magnetoinductance 
L(f) can be found by applying Kirchhoff's laws to the 
inductor network {Lij}, it is evident that, at any frus- 
tration /, L(f) will be completely determined once the 
distribution of {9ij} is known. 

The ground state of a uniformly frustrated array of pe- 
riodically repeated gaskets of order N [see Fig. 1(a)] can 
be constructed by a simple juxtaposition after finding the 
ground state of an isolated iVth-order gasket only if the 
constraints of fluxoid quantization imposed on the largest 
triangular loops (located between the iVth-order gaskets) 
are automatically fulfilled. Recalling the definition of / 
given above, it can be showrJ13 that this condition is sat- 
isfied only for a particular set of frustrations given by: 

f=— (2) 

where M is an integer corresponding to the total number 
of vortices in the rhombohedral unit cell of the SG array 
composed, as shown in Fig. 1(b), by a single iVth-order 
gasket and the adjacent "empty" triangular loop. Thus, 
in order to calculate the magnetoinductance of the sys- 
tem at the values of / given by Eq. (2), all we need is to 
determine the ground-state distribution of in one 

of its constituent gaskets. 

Assuming that the ground state of a gasket is the one 
having the highest possible symmetry consistent with its 
reflection and third-order rotation symmetriesEll (for wire 
networks with the geometry of a third-order SG this con- 
jecture has been_confirmed by numerical calculationsca) , 
it can be showno that the number of independent bond 
variables {Oij} is equal to (3^ + l)/2. They can be found 
from an equal number of constraints imposed by current 
conservation at the nodes and fluxoid quantization in the 
loops of the gasket. For example, in the second-order gas- 
ket of Fig. 2 the links sharing the same symbol have the 
same values of Oij , thereby showing that there are only 5 
independent bond variables. 

Imposing current conservation at one of the 6 equiva- 
lent nodes (denoted by K in Fig. 2) and fluxoid quan- 
tization in the 4 non-equivalent loops a, &, c, and d 
(with, respectively, vortex occupation numbers m a , mj, 
m c , and m,i ) leads to a system of 5 equations (only 
one of which, describing current conservation, is non- 
linear) which have to be solved numerically under the 
constraint that the distribution of the quantum numbers 
{m a } (a — a, b, c, d) is such that the gasket energy: 




FIG. 2. The links of this second-order gasket sharing the 
same symbol have identical gauge-invariant phase differences 
Oij in the states whose symmetry is consistent with the reflec- 
tion and third-order rotation symmetries of the gasket. The 
5 independent {Oij} follow by imposing current conservation 
at one of the 6 equivalent nodes K and fluxoid quantization 
in the 4 non-equivalent loops a, b, c and d. 

£ = J^(l-cos%) (3) 

<ij> 

is minimized. Quite remarkably, for a given frustration 
satisfying Eq. (2), the distribution of {m a } correspond- 
ing to the lowest energy turns out to be identical to that 
emerging from a calculation based on junctions with a 
linear current-phase relation, for which a fully analytical 
treatment is possiblc£irE3. The result is illustrated in Fig. 
3, where the ground-state vortex configurations for the 
rhombohedral unit cell of a regular array of second-order 
gaskets are shown for M in the range [0,40]. 

Inspection of Fig. 3 reveals characteristic features, 
which are valid for arbitrary gasket order. One first 
observes that, with increasing frustration, vortex-jwcle- 
ation spreads from the largest to smallest loopstiE3, a 
property reflecting the hierarchical character of the en- 
ergy needed to create a vortex excitationlI90 Next, one 
notices that the vortices penetrate the gaskets only for 
M > 1, thereby implying that a SG array will exhibit 
fractal behavior only for / > 1/(2 x 4^). One fur- 
ther recognizes that, since for M = 1 the rhombohedral 
unit cell contains just one single vortex sitting in the 
largest triangular loop, the ground state of the array at 
/ = 1/(2 x 4 N ) corresponds to a triangular lattice of vor- 
tices with a nearest-neighbor distance equal to the gasket 
size Ljy. Recalling that Ln plays the role of £ p , one ex- 
pects that for / < 1/(2 x 4^) the system will behave like 
a regular Josephson junction array with conventional eu- 
clidian geometry. Thus, we identify 

fcN = J^N ( 4 ) 

as the frustration at which the EF crossover occurs. The 
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FIG. 3. Ground-state vortex configurations at multiples of 
/ = 1/32 in the range [0,40] for the rhombohedral unit cell 
of a uniformly frustrated periodic array of (undistorted) sec- 
ond-order gaskets. P denotes the number of vortices in the 
largest triangular loop. Vortex configurations with circled 
vortex occupation number P are no longer ground states of 
the weakly distorted gasket of Fig. 6(a) (see Fig. 7 for a 
comparison) . 

nature of the ground states of the SG array in the euclid- 
ian regime (0 < / < /cjy) is discussed in Sec. V. 

Having shown how the structure of the ground state 
can be determined for the particular set of frustrations 
(2), we can now proceed with the calculation of the sheet 
magnetoinductance L(f) of the SG array, the quantity 
measured in our experiments. We first notice that, for 
this particular set of frustrations, L(f) is proportional to 
the magnetoinductance of a constituent gasket. There- 
fore, if we calculate the inductance of a single gasket 
and normalize it to its value at / = to eliminate the 
trivial dependence on the gasket size, we obtain a result 
also expressing the normalized sheet magnetoinductance 
L(f)/L(0) of the composite periodic system. In Ref. [l2| 
we have pointed out that for a given distribution of } 
the calculation of the inductance of a single gasket can 
be performed by successive application of the triangle- 
star transformation well known in the theory of electric 
networkstil. The result of such a calculation for a regu- 
lar array of second-order gaskets is shown in Fig. 4 for 
multiples of / = 1/32 in the range [—11,43]. In order 
to compare this calculation with the experimental data 
presented in the following subsection, we plot the nor- 
malized inverse mag netoinductance L~ 1 (/)/L~ 1 (0). No- 
tice that, as expected for a uniformly frustrated Joseph- 
son junction array, L~ 1 (f) is symmetric with respect to 
/ = 1/2 and periodic with period / = 1. We also recall 
that, although hardly visible in the linear plot of Fig. 4, 
in the fractal regime the frustration-induced correction 
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FIG. 4. Normalized inverse magnetoinductance at multi- 
ples of / = 1/32 in the range [—11,43] of a uniformly frus- 
trated regular array of (undistorted) second-order gaskets. 
Notice the symmetry with respect to / = 1/2 and the pe- 
riodicity of period f — 1. The dotted line is simply a guide 
to the eye. 

AL(/) =.L(f) — L(0) to the array inductance is predicted 
to scaleEj, in the limit of small frustrations and low tem- 
peratures, as f v with v = ln(125/33)/ln4 w 0.96. Obvi- 
ously, the power-law behavior of AL(f) should no longer 
persist in the euclidian regime below f C N. 

So far we have considered only the frustrations given 
by Eq. (2), for which the ground state of a regular array 
of gaskets can be regarded as a periodic replication of the 
ground state of a single gasket. In order to determine the 
array ground state at rational frustrations differing from 
those given by Eq. (2), one should consider a "supercell" 
comprising more than one gasket. The analysis of the 
ground states based on such supercells rapidly becomes 
cumbersome and (for / > / c jv) is not pursued in this 
work. However, as supercells imply that superconducting 
phase coherence extends at larger length scales, one can 
expect the corresponding ground states to be more vul- 
nerable to thermal fluctuations [and, consequently, less 
prominent in the fine structure of L(f)] than those at 
/ = Mf C N- At these particular values of frustration 
the vortex configurations (shown in Fig. 3) are strongly 
pinned by tho-hicrarchical potential landscape provided 
by the gasketscJ, thereby making these ground states par- 
ticularly robust against thermal fluctuations. 

B. Comparison with experiment and effects of 
hidden incommensurability 

Focusing on the fine structure of L(f) we now com- 
pare the theoretical predictions of Sec. Ill A with 
high-resolution magnetoinductance measurements per- 
formed on the array of second-order gaskets. Fig. 
5(a) shows the normalized inverse sheet magnetoinduc- 



tance L _1 (/)/L _1 (0) of the array of second-order gaskets 
(measured at 1 kHz) at three different (reduced) temper- 
atures. We first observe that, although the overall shape 
of the magnetoinductance curves looks roughly similar, 
the fine structure becomes richer and much sharper with 
increasing temperature, thereby revealing very clearly 
almost all the "superfluid" peaks corresponding to the 
states with a unit cell not larger than a single gasket and 
belonging to the sequence / = AT/32 given by Eq. (^). 
Notice that, to make the identification of the structures 
easier, the frustration unit on the horizontal axis of Fig. 
5(a) is chosen to be equal to J c n — 1/32. The striking 
evolution of the fine structure with temperature suggests 
that the motion of vortices due to thermal fluctuations 
plays a major role in the dynamic response of these ar- 
rays. We interpret the behavior shown in Fig. 5(a) as 
clear evidence that, at sufficiently high temperatures, su- 
perconducting phase coherence in the neighborhood of 
the ground states at / = M/32, for which the vortex lat- 
tice is pinned, is drastically disrupted by vortex-lattice 
defects, created by excess or missing vortices, moving 
almost freely on the "frozen" vortex background. This 
process dramatically sharpens the fine structure, thereby 
enhancing the amplitude of the oscillations. Similar be- 
havior was also observed in experimentis-.performed on 
wire networks of interconnected gaskets&a and in regu- 
lar triangular Josephsan junction arraysc^l as well as in 
numerical simulationsEj. 

In sharp contrast to the theoretical prediction for 
a uniformly frustrated SG array (see Fig. 4), the 
curves of Fig. 5(a), although still sym- 
metric with respect to / = 0, are no longer periodic 
with period 1, a behavior leading unavoidably to the 
suppression of the symmetry with respect to / = 1/2. 
We attribute these anomalous features to the inhomo- 
geneous frustration resulting from the change in the ef- 
fective areas of different plaquettes caused by the asym- 
metric (with respect to the direction of the links) dia- 
magnetic response of the truncated-star-shaped super- 
conducting (Pb) islands [see Fig. 1(a)]. Because of this 
particular geometrical form, the screening currents flow- 
ing in these grains create a distortion of the current pat- 
terns associated with the individual loops which leads 
to a redistribution of the magnetic field and perturbs 
the self-similarity of the gaskets. In the temperature 
range (5.5-K" < T < 6AK) of the data shown in Fig. 
5(a), the magnetic penetration depth of the Pb islands 
[as estimated from the zero-temperature bulk Pb value 
(A(0) ~ 40nm) and the (proximity-effect reduced) transi- 
tion temperature (T c = 6.9K) of the Pb grains] is at least 
25 times smaller than their smallest planar geometrical 
dimension, which corresponds to the width of the junc- 
tions (~ 2fj,m). Thus, one expects the distortion of the 
current patterns and, consequently, the nonuniformity of 
the frustration to have a considerable effect on L~ 1 (f). 

The origin of the phenomenon being intimately related 
to the geometry of the superconducting islands rather 
than to the physical properties of the junctions, it seems 
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FIG. 5. (a) Normalized inverse magnetoinductance of the 
periodic array of second-order gaskets shown in Fig. 1(a) 
measured at 1 kHz at three different temperatures as a func- 
tion of frustration. For clarity, the curves for r = 0.13 and 
r = 0.05 are shifted upwards by, respectively, 0.25 and 0.50. 
Notice the absence of periodicity of period 1 and of the re- 
lated symmetry with respect to / = 1/2. (b) Solid circles: 
normalized inverse magnetoinductance calculated for an ar- 
ray of second-order gaskets deformed as shown in Fig. 6 with 
e — 7.4%. This curve should be compared with the data at 
the lowest temperature (r = 0.05) in (a). For comparison, 
the calculation of Fig. 4 for an undistorted SG array is also 
shown (open circles, the curve is shifted upwards by 0.4 for 
clarity). The dotted lines in (b) are simply guides to the eye. 

plausible to describe the effect of the inhomogeneous frus- 
tration by changing merely the effective area of the dif- 
ferent plaquettes according to a prescribed rule. This is 
illustrated in Fig. 6(a) for a second-order gasket, whose 
distortion is modeled by shifting the vertices of the tri- 
angular loops toward the " centers of mass" of the corre- 
sponding superconducting islands, as shown in Fig. 6(b). 
Since the distribution of the screening currents in the 
islands is shifted in the same direction, this deforma- 
tion scheme appears to be a reasonable approach offering, 
above all, the advantage of a simple description in terms 
of a single paramater, the (small) areal change e (ex- 
pressed in units of the area of an elementary triangular 
cell) defined in Fig. 6(a). Notice, incidentally, that this 
procedure does not alter the frustrations [Eq. (2)] for 
which the ground state of the array can be constructed 
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FIG. 6. (a) The deformation scheme (here illustrated for a 
second-order gasket) introduced to describe the effect of the 
nonuniform frustration resulting from the distortion of the 
screening current pattern created by the asymmetric shape of 
the superconducting islands. The vertices of the triangular 
loops are shifted towards the "centers of mass" of the corre- 
sponding superconducting grains as shown in (b). The areal 
changes of the loops a, b, c and d are expressed in terms of 
the deformation parameter e (measured in units of the area 
of an elementary triangular cell) defined in (b). 

by a simple juxtaposition of independent gaskets. 

As the resulting deformation preserves the reflection 
and third-order rotation symmetries of the gasket [see 
Fig. 6(a)], the determination of the vortex and {9ij} 
configurations in the highly symmetric ground states at 
/ = Mf c N can be carried out by following again the pro- 
cedure described in Sec. Ill A. The only modification 
appears in the formulation of the fluxoid-quantization 
constraints for the deformed loops, whose areas (and, 
therefore, the associated magnetic fluxes) change. For 
example, for the second-order gasket of Fig. 6(a) the ar- 
eas of the deformed plaquettes a, b, c, and d are equal, 
in an approximation linear in e, to 1 — e, 1 — (3/2)e, 
1 — 3e, and 4 + 3e, respectively. An analogous problem, 
however with a quite different deformation scheme, was 
considered by Ceccatto et al£3 for a system with a linear 
current-phase relation. 

Ground-state vortex configurations for the rhombohe- 
dral unit cell of an array of second-order gaskets with 
a deformation parameter e = 7.4% (the reason for the 
choice of this value is explained below) are shown in Fig. 
7 for M in the range [0,40]. A comparison with Fig. 3 
reveals that some of the ground states of the undistorted 
array have been replaced by new ones. More precisely, 
while for the regular SG lattice the symmetry with re- 
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FIG. 7. Ground-state vortex configurations at multiples of 
/ = 1/32 in the range [0, 40] for the rhombohedral unit cell of 
a periodic array of second-order gaskets deformed as shown in 
Fig. 6(a) with e = 7.4%. P denotes the number of vortices in 
the largest triangular loop. The comparison with the ground 
states of an array of undistorted gaskets of the same order 
(Fig. 3) reveals that the vortex configurations with circled 
vortex occupation number P are new ground states of the 
system. 

spect to / = 1/2 implies that the ground state for (1 — /) 
may be conceived as resulting from the superposition of 
the ground states for / = 1 and — /, for the distorted 
system this property no longer holds. The absence of pe- 
riodicity (corresponding to the / <-» 1 + / symmetry) and 
of the related / <-> 1 — / symmetry are clearly reflected 
in the normalized inverse magnetoinductance shown in 
Fig. 5(b). This L _1 (/)/-L _1 (0) curve was calculated for 
an array of second-order gaskets using the method out- 
lined in Sec. Ill A and was fitted to the low-temperature 
(t = 0.05) data of Fig. 5(a) using e as an adjustable 
parameter. The best fit was obtained for e — 7.4%. In 
this connection, it should be noticed that, although weak, 
thermal fluctuations still affect the data at r = 0.05 and 
tend to enhance the amplitude of the oscillations with 
respect to that predicted by our calculation, which ne- 
glects thermal fluctuations. Nevertheless, the agreement 
is quite remarkable, especially if one considers that it in- 
volves only one fitting parameter. Moreover, the vertex 
displacement [denoted by the arrow in Fig. 6(b)] defining 
e turns out to be ~ 70% of the distance to the "center 
of mass" of the truncated-star-shaped superconducting 
islands, thereby demonstrating the basic validity of our 
interpretation. 

Similar behavior was observed in the magnetoinduc- 
tance of the array of fourth-order gaskets, whose fine 
structure was found to be richer than that of the array 
of second-order gaskets, as demonstrated by the incip- 
ient splitting of some of the structures at / = M/32. 
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However, because of the strong overlap resulting from 
the nonvanishing width of the superfluid peaks, only a 
fraction of the states at multiples of / = 1/512 could be 
resolved and unambiguously identified in a plot at (rela- 
tively) low frustration resolution like that of Fig. 5 (for 
high-resolution data at very small /, see Fig. 9 in Sec. 
V). 

Before closing this section, we would like to notice that 
when e is a rational number, the periodicity of L~ 1 (f) is 
restored, however with a period larger than 1. 



IV. THE CROSSOVER BETWEEN THE TWO 
REGIMES 

Below f cN [see Eq. (§] the SG arrays are expected 
to be in the euclidian (or homogeneous) regime. In or- 
der to provide preliminary evidence for the existence of 
the EF crossover at f C N, in Fig. 8 we compare, in a 
log-log plot, the quantity AL- 1 (f)/L- 1 {0) = [L _1 (0) - 
L _1 (/)]/L _1 (0), which measures the relative change in 
superfluid density caused by frustration, for the two SG 
arrays studied in this work. Both curves were taken at 
160 Hz and at temperatures such that the structures cor- 
responding to the ordered states are emphasized by ther- 
mal fluctuations. With decreasing frustration the data 
for the sample of fourth-order gaskets exhibit, down to 
/c4 = 1/512, clear fractal features: specifically, four self- 
similar stages (the number of stages being consistent with 
the order of the gaskets) reflecting the dilational symme- 
try of the gaskets and an overall scaling with / which, 
in spite of the indisputable evidence for fluctuation ef- 
fects, follows the asymptotic prediction AL~ 1 (f) oc 
For comparison, the result of a calculation for an undis- 
torted infinite gasket based on the methods discussed in 
Sec. Ill A is also shown in Fig. 8. Below j c \ the data 
tend to flatten out, thereby signaling a possible change 
of regime. However, considering the fact that this change 
sets in almost at the limit of our inductance resolution 
(AL/L ~ 1%), we refrain from drawing a firm conclusion 
with regard to the existence of a dimensional crossover 
in the array of fourth-order gaskets. 

On this subject, the data for the array of second-order 
gaskets convey a much stronger message. In Sec. Ill B 
the magnetoinductance of this sample was shown to obey 
the predictions for the fractal regime whose signatures, 
although less pronounced than in the N — 4 case, can 
also be identified, above f C 2 = 1/32, in the log-log plot of 
Fig. 8 [notice that the self-similar and scaling properties 
of AL~ X are expected to become clearly manifest only-in 
the asymptotic limit (/ — > 0) of higher-order gasketalj]. 
Below f C 2, however, the behavior of AL _1 (/)/L~ 1 (0) 
drastically changes. Besides the loss of self-similarity and 
scaling, the data reveal, by closer inspection, the presence 
of structures [the "dips" in AL~ 1 (f)] corresponding to 
new commensurate states, which can not be ascribed to 
the fractal regime. 
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FIG. 8. Log-log plot of the frustration-induced relative 
change of the inverse sheet inductance of the arrays of sec- 
ond-order (N — 2) and fourth-order (N = 4) gaskets mea- 
sured at 160 Hz. The dotted curve is the theoretical predic- 
tion for an undistorted gasket of infinite order (see Ref. 12). 
The result [the same as in Fig. 5(b)] of the calculation for an 
array of distorted second-order gaskets is shown (open trian- 
gles) at multiples of / = 1/32 in the range [1, 16]. For clarity, 
the theoretical curves are shifted downwards by a quarter of a 
decade. / c jv is the frustration at which the crossover from the 
euclidian (/ < / c jv) to the fractal (/ > / c jv) regime occurs. 



To strengthen the evidence for a regime crossover, in 
Fig. 9 we present the results of sheet inductance mea- 
surements performed at high frustration resolution in the 
range |/| < f c %. Once again, to promote structures corre- 
sponding to ordered states we tuned the effect of thermal 
fluctuations by increasing the temperature and reducing 
the measuring frequency (10 Hz) in an appropriate way. 
Moreover, in Fig. 9 we also include the dissipative com- 
ponent R(f), whose remarkably well resolved fine struc- 
ture provides additional evidence for the EF crossover. 
Notice that maxima in the (normalized) superfluid den- 
sity, L~ 1 (f)/L~ 1 (0), correspond to minima in R(f), as it 
should be. For the array of fourth-order gaskets both the 
superfluid and the dissipative components display clear 
fractal features with marked structures at multiples of 
fd = 1/512 corresponding to the frustration unit on the 
/-axis. In sharp contrast with this behavior, the dy- 
namic response of the array of second-order gaskets is, 
within experimental accuracy, symmetric with respect 
to / = (l/2)/ C 2 and contains well resolved structures 
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FIG. 9. Normalized inverse magnetoinductance and mag- 
neto-resistance of the periodic arrays of second-order (N = 2) 
and fourth-order (N = 4) gaskets measured at 10 Hz in the 
frustration range |/| < f C 2 = 1/32 = 0.03125. Structures 
in the data for N — 2 are labeled in terms of fu = ///c2, 
the frustration parameter referred to the hexagonal unit cell 
of the honeycomb lattice shown in Fig. 10. Structures at 
multiples of f C 4 — 1/512 demonstrate the fractal character of 
the response for N = 4, whereas the prominent structures at 
\fn\ = 1/3 and the symmetry with respect to fu = 1/2 of 
the data for N = 2 are signatures of the euclidian regime. 

at |/| = (l/4)/ c2 , |/| = (l/3)/ c2 and |/| = (l/2)/ c2 . 
In particular, the presence of prominent structures at 
|/| = (l/3)/ C 2 provides compelling evidence for the exis- 
tence of a new regime, since no commensurate state cor- 
responding to this frustration can ever emerge from the 
characteristic sequence / = M J c m of the fractal regime 
for gaskets of any order. We interpret these features as an 
unambiguous signature of the euclidian regime, in which 
the fine structure of the array's impedance reflects, as 
shown in the following section, the existence of vortex 
ground states with a unit cell larger than a second-order 
gasket. 

V. THE EUCLIDIAN REGIME 

In order to understand how our arrays of gaskets be- 
have in the euclidian regime, we first recall that the en- 
ergy of a vortex in a Sierpinski gasket decreases with 
increasing, size of the cell in which the vortex core is 
localizedlld'til. Accordingly, a single vortex in a lattice 
of periodically repeated SGs can be considered as inter- 
acting with an external potential whose minima coincide 
with the centers of the largest triangular cells located 
between the Ath-order gaskets and, therefore, form a 




FIG. 10. Successive applications of the triangle-star trans- 
formation allow to replace each gasket by an elementary 
"star", thereby turning a lattice of periodically repeated gas- 
kets into a honeycomb lattice. 

triangular lattice. For / < / c/ y the number of vortices in 
the system is smaller than the number of the largest tri- 
angular cells, so that it is energetically favorable for the 
vortices to occupy only these largest cells, the concentra- 
tion of the occupied ones being equal to the ratio /// c jv- 
Therefore, one can expect that, for / < / c jv, the behavior 
of an array of periodically repeated SGs resembles that of 
a uniformly frustrated XY-mode\ on a honeycomb lattice 
(with a reduced value of frustration fu = f / f C N < 1), 
whose ground states can also be thought of as formed by 
vortices occupying the sites of a triangular lattice with 
the same concentration fn- 

Another approach leading to the same conclusion relies 
on the iterative procedure, based on successive applica- 
tions of the triangle-star transformations, developed in 
Ref. [l^. At low frustrations, such that there are no vor- 
tices in the gaskets, this method allows to replace each 
gasket by an elementary "star" consisting of three links 
with a modified interaction. This transforms a lattice of 
periodically repeated SGs into a honeycomb lattice (see 
Fig. 10). The difference with respect to the conventional 
XY-mode\ on a honeycomb lattice is that in the iteration 
process the interaction becomes almost harmonir£3, the 
anharmonic corrections becoming smaller and smaller as 
the number of iterations, which is set by the gasket order 
A, increases. Moreover, for / ^ the resulting-pffective 
interaction is no longer an even function of 0^113. 

After establishing that the ground states in the eu- 
clidian regime (/ < /cjv) arc formed by vortices occu- 
pying the sites of a triangular lattice with concentration 
Jh = f '/ 'fcN < 1, it is natural to expect that the states 
which are particularly robust against thermal fluctua- 
tions correspond to values of fu allowing the formation 
of an undistorted triangular vortex lattice (analogous to 
the Abrikosov lattice in bulk superconductors) commen- 
surate with the underlying lattice provided by the array. 
The energy required to create a defect in these highly 
symmetric states is larger than for frustrations requir- 
ing the vortex lattice to be distorted. Therefore, these 
Abrikosov-like states will be less vulnerable to thermal 
fluctuations. For 1/2 < /# < 1 almost equally sta- 
ble states can be constructed when the vacancies in the 
densely packed triangular vortex lattice corresponding to 
fu = 1 also form an undistorted triangular lattice. This 
leads to the symmetry fn <-> 1 — /jj, although it is not 







(a) : f H = 1/3 (b):f H =l/4 




(c):f H =l/2 (d):fg=l/2 (e):f H =l/2 




FIG. 11. (a) and (b) show ground state Abrikosov-like tri- 
angular vortex configurations in a honeycomb lattice at, re- 
spectively, fn = 1/3 and fn = 1/4. (c), (d), and (e) are 
examples of ground-state vortex configurations with the same 
energy in a honeycomb lattice at fn = 1/2. (c): regular ID 
superlattice structure ("striped" phase); (d): a zero-energy 
domain wall separating two "striped" states of the type shown 
in (c); (e): the state obtained from (c) by introducing the 
largest possible number of zero-energy domain walls. 

rigorous. 

It is readily seen that triangular vortex lattices com- 
mensurate with the underlying lattice can be constructed 
for fn = 1/(to 2 + mn + n 2 ), where m and n are integers 
(to > 1, < n < to). In particular, m = n = 1 gives 
fli = 1/3, m = 2 and n — give fjj = 1/4, m = 2 
and n = 1 give fn — 1/7, etc. The states for fn = 1/3 
and Jh — 1/4 are shown, respectively, in Figs. 11(a) and 
11(b). Within this family, the most "dense" state, i.e. 
the state for fn = 1/3, can be expected to be, on ac- 
count of its stronger vortex-vortex interaction, the most 
stable one, followed, in order of decreasing stability, by 
those for Jh = 1/4 and fn = 1/7. This is precisely 
what we observe in the magnetoimpedance data of Fig. 
9 for the array of second-order gaskets where, focusing 
on the interval \,Jh\ < 1/2), we find that the most promi- 
nent structures, signaling a very stable commensurate 
state, appear at \Jh\ = 1/3. Weaker structures corre- 
sponding to the next-stable triangular vortex lattice are 
also well resolved at \fn\ — 1/4, whereas structures at 
\Jh\ = 1/7 are barely visible and present only in the 
inverse magnetoinductance data. It is worth mentioning 
that the observation of the particularly prominent vortex 
state at | /«-!_=. 1/3 is entirely consistent with theoretical 
predictionsLsia as well as with results of Monte Carlo 
simulationsli3 for the frustrated XY model on the hon- 
eycomb lattice. 

The magnetoimpedance data of Fig. 9 show weak 
structures also at \fn\ = 1/2. The ground states of the 
ATF-model on a honeycomb lattice at this particular frus- 
tration and the corresponding vortex configurations were 



studied in Ref. [L4| and shown to possess a so-called ac- 
cidental (i.e., not related to the symmetry) degeneracy, 
which can be discussed in terms of zero-energy domain 
walls. Fig. 11(c) shows an example of a ground state 
for fn = 1/2 characterized by a regular ID superlat- 
tice structure. In this state all the variables take 
the values 6>y = 0,±7r/4. Redistributing the same set of 
variables in a different way among the bonds of the lattice 
allows to transform the state of Fig. 11(c) into another 
one with the same energy, shown in Fig. 11(d), where 
a domain wall separates two " striped" states of the type 
shown in Fig. 11(c). An infinite family of states with the 
same energy can be constructed by creating sequences 
of such zero-energy domain walls parallel to each other. 
Fig. 11(e) shows another example of a periodic ground 
state, which can be obtained by introducing the largest 
possible number of domain walls into the state shown in 
Fig. 11(c). 

It is knowntil that in systems allowing the formation of 
zero-energy domain walls the stabilization of long-range 
order at nonzero temperatures is achieved if the acci- 
dental degeneracy is removed by thermal fluctuations. 
This mechanism (the so-called "order by disorder"E3) 
is relatively weak and, therefore, the ordering in frus- 
trated XF-models with accidental degeneracy should be 
less stable (and destroyed at lower temperatures) than 
in other XY models. This effect should be even more 
pronounced for the fully frustrated (i.e. for fn = 1/2) 
model on a honeycomb lattice, which, in contrast ±o other 
2D AF-models with accidental degeneracyl^JcJTJ, does 
not display any difference in the spin-wave free energy, 
calculated in the harmonic approximation, of different 
periodic ground stateai-J. The Monte-Carlo simulations 
of Shih and Stroudlia have indeed demonstrated that the 
phase transition of the frustrated Xy-model on a hon- 
eycomb lattice at fn = 1/2 takes place at a much lower 
temperature than at fu = 1/3 or fn = 1/4, where the 
accidental degeneracy is absent. 

Returning to the system of periodically repeated SGs 
at / = /cw/2 (Jh = 1/2), it is also possible to start 
the search for its ground state by finding the structure of 
the state which minimizes the energy (||) for the periodic 
vortex configuration shown in Fig. 11(c). Obviously, this 
state is periodic, and its unit cell comprises two gaskets. 
By redistributing the same set of {9ij} among the bonds 
of the lattice in a different way, one can construct the 
state with the same energy corresponding to the vortex 
configuration shown in Fig. 11(d). However, in contrast 
to the fully frustrated XF-model on a honeycomb lattice, 
the state obtained in this way will not be an extremum 
of the Hamiltonian and, therefore, a slight readjustment 
of {Oij} can further decrease its energy. This means that 
in a system of periodically repeated SGs the domain wall 
shown in Fig. 11(d) has a negative energy. Accordingly, 
the state with the lowest energy corresponds to the pe- 
riodic vortex configuration shown in Fig. 11(e), which 
is characterized by the highest possible density of such 
domain walls. The unit cell of this state comprises four 
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gaskets. In the fully frustrated AY-model on a honey- 
comb lattice a ground state with the same structure is 
selected if one takes into account the interaction, of ar- 
bitrary sign, with the second-nearest neighborsEU. 

Thus, the phase transition taking place, with decreas- 
ing temperature, at fjj = 1/2 should be related to the 
appearance of long-range order corresponding to the vor- 
tex configuration shown in Fig. 11(e). The selection of 
this state relies on a weak mechanism, which loses its 
efficiency with increasing N, since the effective interac- 
tion becomes almost harmonic under decimation. One 
can therefore expect the ordered phase at fn = 1/2 to 
be again rather vulnerable to thermal fluctuations. This 
explains why tlm-structures at \Jh\ = 1/2 in Fig. 9 
are much weakera than those at = 1/3. Their 

strength is at most comparable to that of the struc- 
tures corresponding to the triangular vortex lattice at 
fu = 1/4, whose weaker vortex- vortex interaction makes 
the ordering less-, robust than at /# = 1/3. Numer- 
ical simulationsEa show that in the conventional XY 
model on a honeycomb lattice the ordered phases appear, 
with decreasing temperature, in the same order: first at 
fu = 1/3, then at Jh = 1/4, and only further down at 
f H = 1/2. 

VI. CONCLUSION 

Our magnetoinductance measurements on Josephson 
junction arrays of periodically repeated Sierpinski gas- 
kets have clearly demonstrated the existence of two 
regimes. In agreement with theoretical analysis, in one 
of them (the fractal regime) the peaks observed in the 
inverse sheet magnetoinductance, reflecting those states 
which are the most stable against thermal fluctuations, 
are equally spaced. Neighboring states in this sequence 
differ from each other by the penetration of an additional 
vortex into each unit cell of the array. In the other, eu- 
clidian, regime the sequence of the observed stable states 
corresponds to periodic lattices of vortices occupying the 
largest triangular cells of the array. 

The agreement with theory is achieved not only for the 
positions of the different peaks, but also for their relative 
strengths. In the euclidian regime the structure of the 
ordered states is analogous to that on a honeycomb lat- 
tice with reduced frustration and the relative stability of 
the different states can be understood in terms of vortex 
lattice disordering. However, the amplitudes of the peaks 
observed in the fractal regime can be quantitatively ex- 
plained only if the redistribution of the magnetic field in 
the array due to the asymmetric shape of the supercon- 
ducting islands is taken into account. In this connection, 
we would like to mention the recent work by Park and 
HuseE], who compared the energies of different states in 
a wire network with a kagome lattice geometry at full 
frustration. These authors came to the conclusion that 
the effects related to the finite width of the wires can be 



compensated by bending the wires. This is equivalent to 
our conjecture that the influence of the asymmetry as- 
sociated with the screening currents can be reduced to 
a redistribution of the magnetic field with respect to an 
ideal system. 

The results of our magnitoinductance measurements 
show that this phenomenon is more pronounced at high 
frustration levels, an observation consistent with the 
analysis presented in Sec. Ill B. Indeed, in our model, 
the triangular cells exhibiting the largest relative areal 
changes are the smallest (elementary) ones, which are 
precisely those providing the dominant contribution to 
L~ 1 (f) at high values of /. Data taken at very small 
frustrations, like those shown in Fig. 9, are practically 
unaffected by the nonuniform frustration resulting from 
the asymmetry .of the screening currents. 

Early studieso of a large-order (N — 10) single gas- 
ket of superconducting aluminum wires revealed that 
the period associated with adjacent minima of the 
superconducting-to-normal phase boundary T c (/) was 
larger than that extracted from adjacent maxima. We 
have found a similar perturbation of the periodic field 
dependence in our calculations of L _1 (/) based on the 
model for the magnetic field redistribution proposed in 
Sec. Ill B. The same mechanism may be responsible for 
the anomalous feature observed in Ref. 0. It should be 
noticed, however, that the experiments of Ref. [?] were per- 
formed at temperatures very close to T c (0), where the es- 
timated magnetic penetration depth of the Al wires turns 
out to be comparable to their width. Accordingly, the ef- 
fective gasket distortion resulting from screening effects 
is expected to be weaker than in our arrays. A quan- 
titative verification is therefore needed before drawing a 
conclusion as for the ability of our model to explain the 
anomaly observed in Ref. [7| 

The influence of the incommensurability of different 
lattice cells on the magnetoresistance of Josephson junfe 
tion arrays was discussed by Kosterlitz and Granatoo 
in relation to experiments performed on periodic arrays 
with a complex unit cellEi However, quantitative agree- 
ment between experiment and theory in treating incom- 
mensurability effects in Josephson junction prays has 
been demonstrated only by the present workEEl In our 
system the phenomenon of incommensurability is present 
in an hidden form and its manifestations appear as a con- 
sequence of the asymmetric shape of the superconducting 
islands. 
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